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The branching ratio of the r] decay is an important source of information on the value of the 
quark mass ratio = . Furthermore, isospin breaking effects in the decays K ^3n provide 

information on the pion scattering lengths. The cusp effect in the K 3>n decays is presently 
being analyzed by the NA48 and KTeV experiments. From the theoretical point of view, these 
processes have been studied by different methods. We propose a unified and relativistic treatment 
relying on very general principles, unitarity, analyticity and crossing symmetry, combined with 
chiral counting, in order to construct model-independent representations of the corresponding 
amplitudes that are valid at two loops. A general description of the procedure is given and is 
illustrated in the case of the r\ decay amplitude in the leading order in the isospin breaking. 
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1. Introduction 

During the last few years, the decay processes K ^2>n and T] ^ 3?! have been under intensive 
studies, both from the experimental |[l]]-[||] and from the theoretical points of view. The 

importance of these processes, besides the usual determination of decay rates and energy distribu- 
tions, lies in the possibility of studying isospin breaking effects. The appearance of the cusp effect 
in those processes with two neutral pions in the final state enables quite a simple determination of 
the nn scattering lengths (mainly from — > n^n^n^ decay). The T] decays, which are forbidden 
in the isospin limit, offer a good possibility to determine the isospin breaking parameters like R and 
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We present a method how to obtain a model-independent parametrization of the amplitudes of 
the decays 

(1.1a) 
(1.1b) 
(1.1c) 
(l.ld) 

This parameterization hinges on a two-loop construction of these amplitudes based only on general 
properties like unitarity, analyticity, crossing symmetry, relativistic invariance, and chiral power- 
counting for partial wave amplitudes. Up to two loops it takes the form 

J2^{s,t,u) = ^F[^{s,t,u) + '^{s,t,u)] + 0{p^), (1.2) 

where J^p is an overall normalization, while the polynomial part i^{s,t,u) contains free parame- 
ters describing the energy dependence of the processes in analogy to the Dalitz parameters of the 
traditional PDG parameterization. All the non-analytic part of the amplitude connected with the 
final state nn scattering (as discussed below, we could include also other intermediate states but for 
our purpose the tttt ones are enough) is contained in ^ {s,t,u). It depends on the parameters from 
the polynomial part and on parameters describing the nn scattering, like the scattering lengths, 
which in turn allows for their determination from the experimental study of the cusp effects. 

In the following we shall concentrate on the T] decay, where the important physics appears 
already in the first order in isospin breaking (IB^), mainly because of the less involved analytic 
expressions. Naturally, in this limit we cannot describe the cusp effect. The discussion of it is 



postponed to our forthcoming paper []12p, cf. with our previous proceeding [ ]13P where also a small 
summary of the literature concerning it is given. 

The existing theoretical approaches to computing the T] decay amplitudes can be divided into 
several groups. The first one (mainly and [^) encompasses computation using chiral perturba- 
tion theory (ChPT), where the complete two-loop amplitudes at IB^ are obtained. Nowadays, the 
most advanced approach (mainly in the vicinity of the cusp region) describing also the cusp effect 
and already including also the electromagnetic corrections is a method using the non-relativistic 
effective field theory |j5|]. The third group includes the use of analyticity and unitarity for the con- 
struction of the amplitude. In 1996 there appeared almost at the same time two such methods trying 
to find the fixed point solution of Khuri-Treiman type dispersive relations. Kambor et al. [S] and 
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independently Anisowich and Leutwyler solved numerically this integral equation and then ob- 
tained values of the subtraction constants from the matching with ChPT at specific points. These 
approaches perform a numerical resummation of two-pion rescattering contributions. In contrast, 
our method provides an analytic dispersive representation valid to two loops that follows from 
general properties like unitarity, analyticity and crossing symmetry, combined with chiral count- 
ing. For the sake of completeness, let us add that there exist still other approaches using different 
methods, like the apphcation of the Bethe-Salpeter equation to unitarised ChPT [[To[]. 



2. Reconstruction theorem 



As already sketched in our previous proceeding [ 13], we iteratively construct in parallel the rni 



scattering amplitude and the amplitude of Pn Tin, related to the P ^ 371 by crossing symmetry. 



using the same procedure as in [14, 15]. 



Provided we have the following chiral behaviour of the partial waves of the amplitude 

S^{s,t,u) = 1671 ^(/o(5)+ 3/1(5) COS 0)+i2^^^>2, (2.1) 

Re/f=o,i(^)~ 0(7^2)^ lm/,=o,i(^)~0(/), (2.2) 

Rei2/f>2~0(/), Imi;/f>2~0(/), (2.3) 



we can reconstruct the amplitude in the form ( |l.2[ ), with ,'^{s,t,u) being a third order polynomial 
in the Mandelstam variables having the same s,t,u symmetries as the amplitude £/{s,t,u). The 
unitarity part 

'^{s,t,u) =Ws{s,t,u)+WT{t,s,u) + Wu{u,t,s) (2.4) 

is given in terms of single variable dispersive integrals over the imaginary parts of S and P partial 
waves of all the crossed amplitudes. For instance, the contribution of the s-channel amplitude 
AB CD is 

Ws{s,t,u) = 16? ^dx 3 . ^ . [lm/o(x) + 3Im/i(^) ^'"^ 

+A,s\t-u)t^d. j .J:;^;IL^ (2-5) 

thr X[X S) (^)^^/^ (^) 



and similar for the t- and w-crossed-channel contributions Wt and Wy respectively. The triangle 
function Xab{x) is defined as Xab{x) = {x— {m^ + — (m^ — ms)-^). 

To obtain the imaginary parts entering the above expressions, we use the unitarity relation 
projected on the corresponding partial waves, 

lmfpf(s)=l^^^:i^fr'{s) (fl-\s)Yeis-thr,). (2.6) 

The sum goes over all the possible intermediate states k with symmetry factor Sk (in the case of two- 
particle states k, Sk = 2 for undistinguishable states and Sk = 1 otherwise) and thr^; the threshold 
above which this channel opens. In the low-energy region and up to two loops, k are restricted 
to be pairs of light pseudoscalar mesons. In the decay region these can be further restricted to 
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intermediate KK states only. The contributions from other intermediate states, Uke e.g. Kn, can be 
expanded in powers of the Mandelstam variables and absorbed into the polynomial ^{s,t,u). 

Thanks to this restriction on the relevant intermediate states, we can proceed iteratively as is 
shown in the following section on the case of the 77 —>■ n^n^n^ in IB^ But the ai^gumentation goes 
similarly also for the other decays ( |1.1[ ) and the isospin symmetry plays there no important role 
[ jlSl l other than that it simplifies the analytic form of the relations. 



3. Iterative construction of 77 ^ 3;r process 

In the isospin limit the three-pion decay of 77 is forbidden. We shall construct the amplitude of 
this process at first order in isospin breaking, i.e. we shall take the leading order amplitude of this 
decay nonzero, .J/p 7^ 0, and elsewhere in the calculation we will take an exact isospin limit. 

In this limit the different amplitudes are related; the pion-scattering ones (in our sign conven- 
tion) by 

S^c{s,t,u) = —£^x{s,t.,u) — £/x{t,S,u), (3.1) 
£/oo{s,t,u) = —£/xis,t,u) — £/x(,t,S,u) — J2/x{u,t,s) (3.2) 

with the notation 

M)() : TT^Tt" K^K^, S^x ■■ 71° 71+ 71", -K ■ 71+71" ^ TT+TT"; (3.3) 

and the rj decay (or t] tt —> nn) amplitudes by 

J2/Q0{s,t,u) = —£/xiSit,u) — ^/xiti^^u) — £/xi^jt,s), (3.4) 

where 

.s/qo- 1771" ^ 7r°7r°, i4 : r]7r° ^ tt+tt". (3.5) 

It is therefore enough to reconstruct the £/x (reproducing the original Stern et al. computation [|l4|]) 
and s^x amplitudes. 

3.1 Leading order amplitudes 

From ChPT we know that at 0{p^) both amplitudes are represented by first-order polynomials 
in the Mandelstam variables. Their particular choice (connected also with the particular choice 
of the polynomial of the reconstruction theorem) is important since different choices can possibly 
lead to different convergence properties of the chiral expansion and affect the stability of the fit to 



the data. Contrary to | |13| ] we show here the parametrization of s^x using its scattering length and 
effective range parameter. The parametrization of ^ is chosen using the subthreshold parameters. 
The leading-order amplitudes then look like 

£^x{s,t,u) = \6kU + ^{s-AMI)]+0{p^), (3.6) 



Jx{s,t,u)=.^F{AMl+B{s-SQ))+0{p'^), (3.7) 
where = (M^ + 3M^) /3 denotes the center of the Dalitz plot. 
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The 0{p^) chiral perturbation theory result is reproduced by special values of the parameters. 



To get the pion scattering amplitude in terms of the isospin scattering lengths fllq], we can use 

a = 1(^2 _ flO) and b = ^{5a^ - 2a°). The rjn constants would be Jj^AM^ = ^"^^^T^ and 

3.2 First iteration: one-loop expressions 

Using the partial waves computed from these amplitudes inside the reconstruction theorem, 
we obtain the 0{p^) form of the KK scattering amplitude 

£/x{s,t,u) = 1671 {^jcis,t,u) + WAs,t,u) + W+o{t,s,u) + W+Q{u,s,t)) + 0{p^), (3.8) 

where 

W,{s,t,u) = -S7tf{s) (a+^is- 4M2)) (va + ^{s- 20M2)) , (3.9) 

W+o{s,t,u) = - \6nJ{s) (a - + 4^2)) ' + ^ {s - ml) {t - u)^ , (3.10) 

Ji^) = 7^7 f 2 + a(s) In ?^4^) , (3.11) 



16712 V ^ ^ a{s) + lJ' 
a{s) = ^l-^ (3.12) 

and ^^{s,t,u) is a second order polynomial in s, t and u, symmetric mt — u exchange. We can 
write the polynomial in the form 

6 , , A(2) , , , , 

, m) = a + -2 - AMI) + - -^s{s - AMI) " 7^ " ^^1) + «(« " ^^1)] > 

(3.13) 

where 

= - (7a - 16Mi^) - 3^ (69a^ - 456mUj, + 592M1^) (3. 14) 

is the polynomial canceling the real part of the S wave of the unitarity part near the threshold. 
Therefore, the scattering length and the effective range parameter with hats keep their physical 
interpretation up to one loop. Motivated by this interpretation we can resum a part of the terms 
from the two loop chiral order and change also the non-hatted variables into their physical values 
up to one loop. Similarly, we could reconstruct the pion amplitudes up to two-loop order in a 
way where these coefficients keep their physical interpretation up to two loops, just as in the non- 
relativistic approach [^. 

Using the same procedure, we can reconstruct also the rjn^ —>■ k^k^ amplitude, 

s^x{s,t,u) =J^f{^x{s,t,u) + %{s,t,u)+W+Q{t,s,u) + W+Q{u,s,t))+0{p^), (3.15) 
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with 

H',(s,(,«) = -8l/(s)MlW^(7a + 4^(s-6il/5))+S(s-s„)(4a+^(s-12il/;))Y (3.16) 

V^«(V..) = 16./(.) ((^ - AM=) (a- ^<^±|^^) + ^^^^(' - «)) . (3.17) 

=iM2+J(5-5o)+C(5-^o)^ + £>[(?-^o)^ + ("-^o)^]- (3.18) 

Again the hatted coefficients differ from the unhatted ones by the contribution of the 0{p''') order 
(coming from the higher order contributions and the choice of the subtraction scheme). Thus, in 
the fitting of data we can proceed in two ways - either we get the starting value of the unhatted 
coefficients from the fit of the lower order and in the higher order allow only small variations of 
them and thereby actually fit only the polynomial part (eventually together with the nn scattering 
constants); or we can perform a resummation so that we replace the unhatted coefficients by their 
higher order (hatted) values (making an error only of the neglected order - here 0{p^)) and fit just 
this one set of parameters. 



3.3 Second iteration: two-loop result 

Performing the second iteration a few complications occur. In order to get the imaginary parts 
(in fact discontinuities) appearing in the relation (^!^, we need to continue analytically the unitarity 
relation and thereby also the S and P partial wave projections of the NLO amplitudes below their 
physical threshold. This continuation has to be done in correspondence with the right analytic 
structure of the appropriate physical amplitude. This can be obtained by a careful deformation 
of the integration contours in the formulas for the partial wave projections (let us remind that this 
complication is connected with the instability of the particle P and can be understood as an analytic 
continuation in Mp going from the stable values to the instable ones). In the isospin limit this has 
be proven already in JT?]]. Their simple prescription can be used also in the case of P^tt" nn 
beyond the isospin limit - we can avoid intersection of these contours with the cuts attached to 
normal threshold branching points of the computed amplitude and for them there also appears 
no anomalous threshold on the physical sheet. Nevertheless, we believe that the correct analytic 
continuation can be done also for the remaining processes (this will be discussed in our forthcoming 
papers). 

A further problem is related to the fact that the analytic continuation of the discontinuities 
can be complex below the physical threshold (in opposition to the everywhere real discontinuity, 
which equals to double of imaginary part, in the simple case). As a consequence of this and also 
from the contribution of the sunset diagrams (see [^]), the coefficients of the polynomial ^{s,t,u) 
can obtain small imaginary parts. This can be solved again in two ways sketched at the end of 
the previous subsection (by fitting also this small imaginary part of the hatted coefficients; or by 
deviation from the power counting by neglecting this imaginary part). 

In the end we have the two-loop result of the 7771° ^ tt+tt^ amplitude in the form (p.l5[) with 
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an 0{p ) error and 

W,{s,t,u) = ^Y^yj[s) + -^Q^,{s))%{s), (3.19) 

i= 1,2,4,5 ^ ^ J 

(3.20) 

+ F[{t-sof + {u-sof]. 

For the purpose of this proceedings we do not give the explicit form of the polynomials Pi{s) and 
Qi{s) containing the coefficients A, B, C and D and a, b, A^'' and A^^'. The twelve '^i{s) functions 
are two-loop analogues of the one-loop J{s) functions (actually, {s) = J{s)) and of the two- 

loop functions Ki{s) defined in the second paper in [14]. 



4. Conclusion 

We have presented a method that allows to construct a model-independent two-loop parametriza- 
tion of the P ^?>7i decay amplitudes based only on unitarity, analyticity, crossing symmetry, rela- 
tivistic invariance and chiral power-counting for partial wave amplitudes. It contains (polynomial) 
subthreshold parameters describing the energy dependence of these decays and the physical pa- 
rameters of Tin scattering (which can be either obtained from the fit of these processes or one can 
use their values from elsewhere to simplify the fit of the Pn polynomial parameters). We have pre- 
sented the explicit result of this construction for the T] — > 3?! processes in the first order in isospin 
symmetry breaking. 



References 

[1] A. Abele et al, Phys. Lett. B 417 (1998) 197; 

F. Ambrosino et al, J. High Energy Phys. 0805 (2008) 006 [ai-Xiv:0801.2642 [hep-ex]]. 

[2] A. Abele et al., Phys. Lett. B 417 (1998) 193; 
M.N. Achasov etai, JETP Lett. 73 (200 1 ) 45 1 ; 
W. B. Tippens et al., Phys. Rev. Lett. 87 (2001) 192001; 
F. Ambrosino et al., arXiv:0707.4137 [hep-ex]; 

M. Bashkanov et al, Phys. Rev. C 76 (2007) 048201 [arXiv:0708.2014 [nucl-ex]]; 
C. Adolph et al, arXiv:0811.2763 [nucl-ex]; 

S. Prakhov et al, Phys. Rev. C 79 (2009) 035204 [arXiv:0812.1999 [hep-ex]]; 
M. Unverzagt et al., Eur. Phys. J. A 39 (2009) 169 [arXiv:0812.3324 [hep-ex]]. 



7 



Dispersive construction of two-loop P ^3n {P = K,r\) amplitudes 



Martin ZDRAHAL 



[3] J. R. Batley etal. [NA48/2 Collaboration], Phys. Lett. B 633 (2006) 173 [arXiv:hep-ex/0511056]; 

NA48 talks [see http://na48.web.cern.ch/NA48AVelcome/images/talks.htmll; 

E. Abouzaid et at. [KTeV Collaboration], Phys. Rev. D 78 (2008) 032009 LarXiv:0806.3535 [hep-ex]]. 

[4] J. Bijnens and F. Borg, Nucl. Phys. B 697 (2004) 319 [arXiv:hep-ph/0405025]; 
N. Cabibbo, Phys. Rev. Lett. 93 (2004) 121801; 

N. Cabibbo and G. Isidori, J. ffigh Energy Phys. 0503 (2005) 021 [arXiv:hep-ph/0502130]; 
E. Ganniz, J. Prades and L Scimemi, Eur. Phys. J. C 50 (2007) 405 [arXiv:hep-ph/0602023]. 

[5] G. Colangelo, J. Gasser, B. Kubis and A. Rusetsky, Phys. Lett. B 638 (2006) 187 
[arXiv:hep-ph/0604084]; 

M. Bissegger et ai, Phys. Lett. B 659 (2008) 576 [arXiv:07 10.4456 [hep-ph]] and Nucl. Phys. B 806 
(2009) 178 [arXiv:0807.0515 [hep-ph]]. 

[6] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 539. 

[7] J. Bijnens and K. Ghorbani, J. ffigh Energy Phys. 0711 (2007) 030 [arXiv:0709.0230 [hep-ph]]. 

[8] J. Kambor, C. Wiesendanger and D. Wyler, Nucl. Phys. B 465 (1996) 215 [hep-ph/9509374]. 

[9] A. V. Anisovich and H. Leutwyler, Phys. Lett. B 375 (1996) 335 [hep-ph/9601237]. 
[10] B. Borasoy and R. Nissler, Eur. Phys. J. A 26 (2005) 383 [hep-ph/0510384]. 
[11] J. Bijnens, PoS (EFT09) 022 [arXiv:0904.3713 [hep-ph]]. 
[12] K. Kampf, M. Knecht, J. Novotny, M. Zdrahal: in preparation. 

[13] K. Kampf, M. Knecht, J. Novotny and JVI. Zdrahal, Nucl. Phys. Proc. Suppl. 186, 334 (2009) 
[arXiv:0810.1906 [hep-ph]]. 

[14] J. Stern, H. Sazdjian and N. H. Fuchs, Phys. Rev. D 47 (1993) 3814 [arXiv:hep-ph/9301244]; 
IVI. Knecht, B. IVIoussallam, J. Stem and N. H. Fuchs, Nucl. Phys. B 457 (1995) 513. 

[15] M. Zdrahal and J. Novotny, Phys. Rev. D 78, 116016 (2008) [arXiv:0806.4529 [hep-ph]]. 

[16] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

[17] J. B. Bronzan and C. Kacser, Phys. Rev. 132 (1963) 2703; 

V. V. Anisovich and A. A. Anselm, Usp. Fyz. Nauk 88 (1966) 287. 



8 



